Quantitative analysis of directional spontaneous emission spectra from light sources in 

photonic crystals 



Ivan S. NikolaevQ Peter Lodahl, and Willem L. Vos 

Complex Photonic Systems (COPS), Department of Science and Technology, 
and MESA + Institute of Nanotechnology, University of Twente, 
PO Box 217, 7500 AE Enschede, The Netherlands^ 

We have performed angle-resolved measurements of spontaneous-emission spectra from laser dyes 
and quantum dots in opal and inverse opal photonic crystals. Pronounced directional dependencies 
of the emission spectra are observed: angular ranges of strongly reduced emission adjoin with 
angular ranges of enhanced emission. It appears that emission from embedded light sources is 
affected both by the periodicity and by the structural imperfections of the crystals: the photons 
are Bragg diffracted by lattice planes and scattered by unavoidable structural disorder. Using a 
model comprising diffuse light transport and photonic band structure, we quantitatively explain 
the directional emission spectra. This work provides detailed understanding of the transport of 
spontaneously emitted light in real photonic crystals, which is essential in the interpretation of 
quantum-optics in photonic band-gap crystals and for applications wherein directional emission and 
total emission power are controlled. 
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I. INTRODUCTION 

Photonic crystals attract much attention both in 
academia and in industry because they offer exciting 
ways of manipulating photons pj. Periodic variations of 
the refractive index in photonic crystals on a length scale 
of the wavelength of light cause optical Bragg diffraction 
and organize the photon dispersion relation in bands, 
analogous to electron bands in semiconductors 0, El E| 
Frequency windows, called stop bands, appear in which 
there are no modes for certain propagation directions. 
Photonic structures already serve as a base for control- 
ling the propagation of light. Of even greater interest 
are three-dimensional (3D) photonic crystals that sum- 
mon novel opportunities in the case of a photonic band 
gap (PBG) - a frequency range where no modes exist at 
all. The existence of a gap in the density of photonic 
modes leads to novel quantum-optical phenomena such 
as spontaneous emission inhibition and light localization 

[limn. 

Control over the radiative decay rate of spontaneous 
emission is of keen interest for applications, and therefore 
emission properties of sources such as atoms, dyes, and 
quantum dots are intensively investigated. According to 
Fermi's 'Golden Rule', this decay rate is proportional to 
the local radiative density of states (LDOS) that counts 
the number of electromagnetic modes wherein photons 
can be emitted at a specific location of the source. In 
3D photonic crystals pronounced variations of the LDOS 
with frequency are predicted even in the absence of a 
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which give rise to angle-independent varia- 
tions of the total emission. Recently LDOS effects on 
spontaneous emission have been experimentally demon- 
strated: Using inverse opals photonic crystals consid- 
erable variations of the emission rates in large band 
widths were obtained in both continuous- wave (cw) total- 
emission power ex peri ments |9| and in time-resolved life- 
time experiments |l(lj | . While lifetime experiments pro- 
vide a direct measurement of decay rates, it is important 
to quantitatively interpret concomitant emission spectra, 
for instance, to confirm that light sources inside the crys- 
tal are probed. Cw experiments, on the other hand, rely 
on a comparison of angle-integrated spectra with a ho- 
mogeneous medium. In the latter case a complete un- 
derstanding of all angle-dependent effects, that is, Bragg 
diffraction on the propagation of light is crucial. In this 
paper such a quantitative analysis is presented. 

In contrast to decay rates, emission spectra of sources 
embedded in pho tonic crystals are strongly directional 
[ill Ha . Il3 . [Tj, [i^, O ■ Particular frequency ranges of the 
spectra are suppressed in certain directions forming stop 
bands, whose center frequencies and the widths are de- 
scribed by the photon dispersion relation. Besides Bragg 
diffraction, which is an effect of the order of the periodic 
structure, light propagating inside the structure also feels 
disorder: polydispersity, roughness and misarrangements 
of the building blocks [17| • This unavoidable disorder af- 
fects the interference-induced properties of photonic crys- 
tals. Previous work on the effect of disorder on sponta- 
neous emission includes the realization that disorder de- 
termines the depth of stop bands [u3 . Il5| . Furthermore, 
the first observations of enhanced emission in the range 
of first and second-order stop bands were also related to 
disorder- induced redistribution of emitted photons [l8|| . 
While the propagation of light from external light sources 
has been studied in great detail 0, [23, Hl| , a quantita- 
tive explanation of the behavior of light emission from 
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internal sources has lacked sofar. 

Here we present strongly frequency-dependent angular 
distributions of spontaneous emission from a laser dye in 
polystyrene opals and from quantum dots in titania in- 
verse opals in the frequency range around the first-order 
Bragg diffraction (L-gap). We compare the data to a the- 
oretical model that unifies effects of structural disorder 
and photonic crystal properties plj . Angle-dependent 
internal reflection due to the photonic gaps plays a key 
role in our model. The theory quantitatively explains 
both the enhancement and the reduction of light along 
certain propagation directions that were observed exper- 
imentally. The excellent agreement confirms that the 
propagation of light in a photonic crystal is well under- 
stood for frequencies around the L-gap. Furthermore, we 
show that by analyzing the exit emission distributions, 
one can reveal stop bands in the quantum dot spectra. 
Such an analysis should be carefully performed before 
any quantum-optical experiments since it unambiguously 
proves the effect of the photonic crystal on emission. We 
finally discuss the applicability of photonic crystals for 
improvement of the emission efficiency of light sources. 



II. EXPERIMENTAL DETAILS 
A. Samples 

We have studied emission from dyes in polystyrene 
opals and from quantum dots in titania inverse opals. 
The polystyrene opals are fee crystals of close-packed 
polystyrene spheres prepared from a colloidal suspen- 
sion by self-assembly. The titania inverse opals are fee 
structures of close-packed air spheres in a solid matrix 
of TiC>2. Details of the preparation and characterization 
of the opals and inverse opals can be found in Ref . [22( • 
We have studied 4 polystyrene samples with lattice pa- 
rameters a — 178 ± 3 and 365 ± 5 nm and 8 titania 
samples with lattice parameters a — 370, 420, 500, 580, 
and 650 ± 10 nm. All samples have typical dimensions 
of 2x2x0.2 mm 3 and contain high-quality domains with 
diameters larger than 50 fim. These domains have flat 
faces normal to the 111 crystal axis, which is evident 
from SEM images. The other crystalline axes are ran- 
domly oriented. 

The polystyrene opals were doped with the laser dye 
Rhodamine 6G (R6G) by soaking them for 30 minutes in 
a dilute solution (10~ 6 mol/1) of the dye in ethanol. Af- 
terwards the samples were rinsed in ethanol and dried 
to remove the dye from the sample surface. To esti- 
mate an upper bound to the density of the dye, we as- 
sume that the infiltrating solution completely fills the air 
voids in the opals, and that, in the process of drying, 
the dye molecules distribute uniformly on the surfaces 
of the spheres inside the opals. Knowing the volume 
of the air voids in each unit cell of the opals (0.26 a 3 ) 
and the dye concentration, we arrive at no more than 10 
dye molecules per unit-cell inner surface |23j| (for opals 



with lattice parameter a — 365 nm). While this surface 
density increases linearly with the lattice parameter a, 
the average distance between the dye molecules remains 
more than 2 orders of magnitude away from the typical 
intermolecular distances where reabsorption and energy- 
transfer processes could play a role [24[ • Before emission 
experiments, we performed a selective bleaching 251 of 
the dye at the surface of opal photonic crystals in or- 
der to ensure that the emission is recorded only from the 
bulk of the crystal, and not from the crystal surface ^| . 
The whole surface was illuminated for up to 1 hour by an 
intense laser beam at the Bragg angle for the pump fre- 
quency. At this angle the pump intensity decreases expo- 
nentially with depth, which implies that the dye bleaches 
within the first few crystal layers. However, our exper- 
imental results were found to be almost independent of 
the amount of bleaching, thus indicating that emitters in 
the bulk of the photonic crystals provide the dominating 
contribution to the measured emission intensity. 

The titania inverse opals were infiltrated for 24 hours 
with a colloidal suspension of ZnSe-coated CdSe quantum 
dots Ha,|27| in a mixture of 50% chloroform and 50% bu- 
tanol. Afterwards the samples were rinsed in chloroform 
and dried. The concentration of the infiltrating solution 
was 10 -7 mol/1. We used the same estimation in order 
to get the maximum number of quantum dots per unit 
cell in the inverse opals. For samples with the largest 
lattice parameter a — 650 nm, this concentration is 15 
quantum dots per unit cell and is sufficiently low to avoid 
energy-transfer processes and reabsorption. In order to 
minimize oxidation and contamination of the quantum 
dots, the inverse opals were infiltrated with the quantum 
dots in a nitrogen-purged glove box and held in a sealed 
chamber under a 1.7 mbar nitrogen atmosphere during 
the optical measurements. No bleaching at the surface 
layers could be performed for the quantum dots, however 
a detailed analysis of the angular emission data (see be- 
low) demonstrates that light emission from the quantum 
dots from the bulk dominates. 



B. Experimental set-up 

Figure ^ shows the experimental set-up used to mea- 
sure emission from light sources inside photonic crystals. 
The sources inside the crystal are excited by a cw Ar-ion 
laser (A = 497 nm) with the power at the sample around 
10 fiW. At this pump power, we do not observe any ef- 
fects of bleaching of the dye during the emission experi- 
ments. To focus the pump beam on the sample surface 
a fiber-coupled microscope objective is used. The beam 
is focused to a spot of about 30 /an in diameter at an 
incident angle 9 p relative to the surface normal, usually 
9p « 25°. In order to acquire emission spectra as a 
function of the detection angle 6 e relative to the surface 
normal, the sample is mounted on a rotation stage. The 
surface normal is parallel to the 111 reciprocal lattice vec- 
tor Gin. In order to illuminate the same area irrespec- 
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FIG. 1: (color online) Scheme of the experimental set-up. The 
pump beam is focused onto the sample at incident angle 9 P 
by the objective 0\ (f = 7.3 cm, NA = 0.05). Luminescence 
within a cone centered at detection angle 9 e relative to the 
surface normal is collected by the lens L\ (f = 12 cm) and 
imaged on the spectrometer slit by the lens L2 (f = 12 cm). 
A colour filter F prevents scattered pump light from entering 
the prism spectrometer. The angle 6 e is varied by rotating 
the rotation stage, which carries the sample holder, the fiber, 
and the objective 0\, whereas the incident angle 6 V is kept 
fixed. 



tive of 8 e , the fiber-coupled objective {0\) is mounted 
on the same rotation stage as the sample. In this way, 
the angle of incidence 9 p between the pump beam and 
the surface normal is kept constant. The advantage over 
previous experiments, where the sample was rotated with 
respect to both the pump and detection beams, is that 
the absolute intensity of the angle-dependent spectra can 
be reliably compared. The position of the pump spot on 
the sample is monitored with a microscope. The emitted 
light is collected within a cone of 15° full width around 
the angle 8 e , imaged on the slit of the spectrometer with 4 
nm resolution, and detected by the photomultiplier tube 
(PMT). The angle-resolved spectra are usually measured 
at the detection angles from 0° to 75° at intervals of 15°. 
The measured spectra are corrected for the dark count 
of the PMT. The shapes of the spectra are confirmed to 
be independent of the pump intensity, and the emitted 
intensity is linear with the pump power. 



III. DIFFUSE LIGHT TRANSPORT IN 
PHOTONIC CRYSTALS 

A. Escape function 

In real photonic structures, defects in the arrangement 
of the building blocks are always present and cause ran- 
dom multiple scattering of light. This means that all 
light emitted in such photonic structures becomes diffuse 
on length scales equal to the transport mean free path 
I, which is often much smaller than the thickness of the 
sample L. For example, our opals and inverse opals have 
mean free paths of about 15 /jm 28], whereas the thick- 



ness of the samples is about 200 /xm. Thus, even though 
photons generated inside a photonic crystal are diffracted 
by the crystal structure, this effect is smeared out by the 
random multiple scattering while the photons propagate 
through the bulk towards the crystal surface. Only at 
distances from the surface z smaller than I, where the 
photons emanate ballistically towards the crystal-air in- 
terface after a last scattering event, the effect of Bragg 
diffraction is not destroyed by the scattering. Hence the 
diffuse emission acquires a directional dependence only 
when it exits the crystal plf . 

We consider the ratio of the mean free path / to the at- 
tenuation length for Bragg diffraction Lb in order to esti- 
mate the attenuation of emission caused by Bragg diffrac- 
tion, as proposed in Ref . 0] . Since the mean free path I 
is larger than the Bragg attenuation length Lb (typically 
1/ Lb ~ 2 — 5 an attenuation in the stop band equal 

to 1 — Lb /I = 50% to 80% is predicted, which is in agree- 
ment with our observations. As will be discussed later, 
the stop-band attenuations are obtained directly from re- 
flectivity measurements, therefore, the mean free path I 
is not an explicit parameter in our theoretical model. 

In the present work we investigate directional proper- 
ties of light emitted by sources from 3D photonic crys- 
tals and compare to a model of diffuse light transmis- 
sion through opaque media [2^, |3(], extended to pho- 
tonic crystals |21|. Based on the diffusion theory, the 
intensity of light I(oj, fi e ) with frequency w that exits the 
sample at external angles between 9 e = cos (/i e ) and 
cos _1 (/i e + d/i e ) relative to the surface normal is equal to 



I(oj, 



ll e )dn e = It ot (u>)P(cO,fX e )dfJ, e 



(1) 



Here L to t(to) is the total spontaneous emission power that 
is the spectrum of the light sources integrated over the 
exit angles 9 e . For sources with a low quantum efficiency 
or with inhomogeneously broadened spectra, Itot(w) is 
proportional to the LDOS The distribution P(w,/i e ) 
is defined as 



nl(l + R(u) 



1 - R(u) 



+ - f i i )[l-R(to, tM )], (2) 



where n e and rii are average refractive indices out- 
side and inside the sample 32], respectively. /i e and 
fii are related by Snel's law. R(u),iJ,i) is an angle- 
dependent internal-reflection coefficient that yields an 
angle-averaged internal-reflection coefficient R(u>): 



- = 3C 2 H + 2C 1 H 
1 ' 3C 2 (uj) - 2Ci(w) + 2 



C n {uj) 



R 



(3) 
(4) 



From the diffusion theory, R{u>) determines the so-called 
extrapolation length that sets the boundary conditions of 
the diffuse intensity |29ll3*cl l3*H . The normalized function 
P(u>, fJ. e ) describes the distribution of emission intensity 
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over the escape angles and will be called the 'escape func- 
tion'. In absence of reflection effects, the escape distibu- 
tion tends to the well-known Lambertian distribution of 
diffuse surfaces. 

In random media such as powders or macrop- 
orous sponges the internal-reflection coefficient R(u>, /it) 
is barely frequency dependent |33l |. and propagation 
through the interface is well described by Fresnel re- 
flection model assuming an average refractive index. 
The angular dependence of the escape function P(u>. Ue) 
agrees well with experiments on random media 30, 31]. 
For highly dispersive photonic crystals, however, Fresnel 
model cannot describe the internal reflection since light 
escaping from a depth z < I from the crystal surface is 
Bragg attenuated for angles and frequencies inside a stop 
band. We model the strong angle and frequency depen- 
dent internal reflection with photonic band structures. 
At a particular frequency uj* where a stop band is present, 
the internal reflectivity R(oj*,[ii) blocks the emission in 
the directions of the stop band (a range of fi* 's related to 
lu* by the photonic band structure) and therefore reduces 
the escape function P(oj*, ^*), c/. Eq. (J2J. The presence 
of the stop band raises the angle-integrated reflectivity 
R(u)*), which enhances the escape function P(ui*,u. e ) for 
angles outside the stop band. Thus, the escape function 
is strongly non-Lambertian in a photonic crystal, show- 
ing clear suppressions or enhancements. 



B. Internal- reflection coefficient 

In order to model the internal-reflection coefficient 
we have taken into consideration calculated 
photonic band structures. Figure Ela) shows the pho- 
tonic band structure for polystyrene opals calculated 
along the LU and LK lines in reciprocal space in the 
frequency range around the first-order stop gap (L-gap) 
that is due to optical Bragg diffraction by (111) planes 
parallel to the sample surface. The emission spectrum 
of R6G is in the low-frequency limit relative to the stop 
bands of opals with a lattice parameter a = 178 nm (re- 
gion confined by the orange dash-dot-dotted lines) . Con- 
sequently, this sample is effectively homogeneous for the 
emission frequencies, and therefore it can serve as a re- 
liable reference. For these nonphotonic crystals we used 
Fresnel model in order to describe the internal reflection. 
In Figure [3f a) one can also see the frequency gap be- 
tween the two lowest bands (blue dashed curves), which 
obeys Bragg's law (red curve) within the frequency range 
of R6G emission for the opals with a lattice parameter 
a = 365 nm (region confined by the green dash-dotted 
lines). Therefore, the angular dependence of the center 
frequency of the L-gap is modelled with the red curve, 
i.e.: uc(/ii) = ^>c{^i — 1)/Mi- To investigate the fre- 
quency dependence of the reflectivity, we have performed 
reflectivity experiments on the samples using external in- 
cident plane waves, since this technique reveals the center 
frequencies and the widths of stop bands [20I l3^ | . Fig- 



ure [3b) shows a normal-incidence reflectivity spectrum 
measured from an opal with the lattice parameter a = 
365 nm (blue circles). The reflectivity peak is not fitted 
well with a Gaussian (red dashed curve). As an improved 
model, we propose a modified Gaussian (Fig. Efb), black 
curve) : 



Ri(cj,fj,i) = Ai(^)exp 



2(A^ c (/i. i )) 4 



(5) 



where Ai(/ij) is the magnitude of the internal- reflection 
coefficient and Awc(A t i) is the width parameter. This 
peak shape is seen to fit the measurements well for fre- 
quencies > 16000 cm -1 . At frequencies below the stop 
band, i.e. below 16000 cm -1 for these particular samples 
with a = 365 nm, a deviation from the model is observed. 
We attribute this deviation to Fresnel reflection, which 
is important only in the low-frequency limit and there- 
fore is not relevant for the escape function of photonic 
samples. The width of the L-gap AujciHi) hardly varies 
with yu,; within the range of the dye emission, therefore it 
is taken to be constant in our model. The magnitude of 
the internal-reflection coefficient A\ (/i^) decreases with /x.; 
because at larger internal angles 0i — cos _1 (/ii) the path 
length for the light to become Bragg attenuated increases 
with /Zi, and this increases the probability of scattering 
at z < Lb- The value of Ai(/j,i) at Hi = 1 is taken from 
the normal-incidence reflectivity experiments. Thus, we 
have: Ai(fii) = Ai(ui = 1) >m and A\{jii = 1) = 0.7, see 
Fig. lib). 

Emitted light that is scattered within a distance Lb < 
z < I towards the exit interface can also be redirected by 
Bragg diffraction by the sets of (111) planes, which are 
oriented at 9 L = 70.5° to the (111) planes and the surface 
normal. The internal-reflection coefficient i?2(^,/^i) for 
Bragg diffraction by (111) lattice planes is modelled sim- 
ilarly to i?i(w, /ij). Taking into account that we measure 
emission from many randomly-oriented crystal domains 
in azimuthal directions, the reflectivity i?2(^,$i) is av- 
eraged over the azimuthal angles <f> between the LK and 
LU lines in reciprocal space, yielding: 



7T 



7r/3 



Ao 



i exp 



(6) 

The magnitude A 2 (9i,<j>) is modelled as A 2 (6i,(j)) — 
A 2 (70.5°, 0°) • cos(9i - 70.5°) • cos{cf>) with A 2 (70.5°, 0°) = 
0.7. The total internal-reflection coefficient R(uj,fj,i) is 
calculated as a sum of the i?i(w,/Zi) and R2(Lo,cos(9i — 
70.5°)) modified Gaussian peaks. We expect this model 
of the angle- and frequency-dependent internal reflectiv- 
ity to capture the essential frequency dependence of the 
first-order photonic stop bands in polystyrene opals. 

In the case of the titania inverse opals we apply the 
same escape model to explain our experimental data. 
However, in calculating the internal-reflection coefficient 
R(to,lii), Bragg diffraction from other lattice planes 
must also be included. This difference compared to 
polystyrene opals appears since titania inverse opals are 
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External angle 9 e [deg] Frequency [cm 1 ] 

FIG. 2: (color online) (a) Photonic band structure for polystyrene opals (blue dashed curves), center frequency of the stop 
band vs. detection angle 8 e according to Bragg's law (red curve), the frequency range of R6G emission is shown for opals with 
lattice parameters a = 178 and 365 nm (hatched regions between the orange dash-dot-dotted lines and the green dash-dotted 
lines, respectively), (b) Normal-incidence reflectivity as a function of frequency for polystyrene opals with a lattice parameter 
a = 365 nm and refractive-index contrast m = 1.59. Blue circles are measured values, the black curve is a fit with the modified 
Gaussian reflectivity peak (Eq. the red dashed curve is a fit with a Gaussian reflectivity peak. 



more strongly photonic and the measurements were per- 
formed at higher reduced frequencies (a/ A = Loa/2-Kc). 
Moreover, the resulting stop bands occur at lower detec- 
tion angles 9 e in these crystals than in the polystyrene 
opals, as a consequence of the lower average refractive 
index. Therefore Bragg's law is not a sufficient approxi- 
mation and is not used to model the angular dependence 
of the stop bands in the titania inverse opals. Instead, 
the full band structure model is employed, in which we 
take into account multiple Bragg wave coupling 0, |3(| . 
For the inverse opals, this model was already successfully 
tested on diffuse transmission experiments [21| . 



IV. RESULTS AND DISCUSSION 

A. Spontaneous emission of R6G in polystyrene 
opals 

Reflectivity measurements at normal incidence of 
polystyrene opals (Fig. |2{b)) reveal that the relative 
width of the first-order stop band is Auj/lo » 0.075. For 
opals with a lattice parameter a — 365 nm this means a 
stop band in the range 16100 - 17300 cm" 1 for light es- 
caping the crystal normally to the surface. The dye R6G 
emits in the range of 15000 - 20000 cm -1 , and hence we 
expect to observe directional-dependent emission of the 
dye from the opals with a = 365 nm. Figure^a) displays 
the emission spectra at selected detection angles for such 
doped opals. It is clearly seen that the shapes of the spec- 
tra are affected by the photonic crystal. The emission is 
suppressed by the crystal for 9 e = 0° in the spectral 
range from 16000 to 17500 cm -1 . With increasing angle 
9 e the low-frequency parts of the emission recover, and 
the suppressed emission range shifts to higher frequen- 
cies, as expected from Bragg's law for a photonic stop 
band from a single set of lattice planes. In contrast, the 
shape of the spectra from an opal with a lattice parameter 



a = 178 nm remains unchanged (Fig.yjb)). The sample 
is not photonic for the frequency range considered: the 
frequencies of R6G emission lie far below the first-order 
stop band in the opal with this lattice parameter. 

Before studying spontaneous emission from photonic 
samples we have verified the applicability of the above- 
mentioned model of diffuse light transport on the non- 
photonic, reference samples. We use Fresnel reflection to 
describe the angular-dependent internal-reflection coef- 
ficient, taking an average refractive index n av = 1.44, 
which is derived from the polystyrene filling fraction 
tp w 74% in opals and the refractive index of polystyrene 
n = 1.59. We record the intensity at the maximum of 
the emission spectrum as a function of the exit angle 9 e 
relative to the measurement at 9 e = 0°. The relative 
intensity I{9 e ) is compared to the escape function P(9 e ) 
in Figure 01 While the expected decrease with angle is 
observed, it is clear that the calculated intensity differs 
systematically from the measured data. This deviation 
appears to be caused by an angle-dependent detection 
efficiency as a result of an increase with 9 e of the pro- 
jection of the spectrometer slit on the sample. Correct- 
ing the measured intensity I{to, fi e ) for the detection effi- 
ciency D(fj, e ) (see Appendix) yields the corrected inten- 
sity I c (w, [i e ) = I(uj, ^i e )/Z)(/i e ) displayed as red symbols 
in Fig. 2| The agreement between the corrected intensity 
and the calculated escape function is excellent. With the 
proper account of the detection efficiency, the angular 
distribution of emission escaping the reference samples is 
thus fully understood. In all experimental data presented 
in the remainder of this paper the detection efficiency has 
been included. 

In the case of the photonic samples the exit distribu- 
tion of emission strongly depends on the frequency u) as 
mentioned above: P(uj,9 £ ) — I c (LU,9 e )/I to t{uj)- The to- 
tal emission spectrum I to t(u) is determined by discretely 
summing the angle-resolved spectra I c (u>,9 e ) weighted 
by 2-7T sm(9 e )d9 e to approximate the integration over 27T 
solid angle. The spectra from Figure 01a) divided by the 
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FIG. 3: (color online) Emission spectra of R6G in polystyrene opals with lattice parameters a — 365 nm (a) and a — 178 nm 
(b). The black curves are obtained at 8 e = 0° , the red dashed curves at 6 e = 30°, the green dotted curves at <9 e = 37.5°, and 
the blue dash-dotted curves at 8 e = 60°. The magenta short-dashed curve in (a) indicates the total emission spectrum 7 tot (oj). 



total emission spectrum Itot(u>) are plotted in Figure 
(symbols) together with the calculated escape function 
P(w, e ) (curves). We observe a good agreement between 
our experiment and theory. The escape function hardly 
varies with frequency in the low-frequency region < 15600 
cm , while it still depends on the detection direction. In 
contrast, at higher frequencies strong variations are seen 
compared to the low-frequency range. At the exit angle 
6 e = 0°, the escape function is significantly reduced in 
the spectral range from 16000 to 17500 cm -1 by the stop 
band centered at u> = 16700 cm -1 due to internal Bragg 
diffraction, which is described by the term (1 — R(u>, fii)) 
in Eq. J2J. The change of the center frequency as well 
as the decrease in the attenuation of emission inside the 
stop band with increasing exit angle 9 e are well described 
in our model by the frequency and angular dependent 
internal-reflection coefficient R(uj,Hi). At 6 e = 60°, the 
stop band has moved out of the spectral range of R6G. 

Figure El also shows a peculiar feature: the frequency 
ranges where the emission is inhibited along certain direc- 
tions, adjoin with the ranges where emission is increased 
along the same directions. Such an increase appears at 
the blue side of the stop band at 8 e = 0° and 30°, and 
at the red side of the stop band at 6 e — 45° and 60°. 
This enhanced escape probability in the frequency range 
16000 - 19500 cm -1 along directions that do not coin- 
cide with the stop band is reflected in our model for 
P(w,/i e ) (Eq. @) as an increase of the angle-averaged 
internal-reflection coefficient R(ui). The good agreement 
between experiments and theory confirms a qualitative 
attribution of such enhancements to diffuse escape ef- 
fects [l8| . Moreover it unambiguously demonstrates that 
our experimental observation of the emission enhance- 
ment is not due to Bragg standing wave effects proposed 
in Ref. j3^, but is related to diffusion of light. In 
Figure we compare the experimentally determined in- 
tensity distributions (symbols) for several fixed frequen- 
cies with the calculated ones (solid curves) as functions 
of the exit angle 9 e . The experimental values of P(0 e ) 
were obtained by dividing the emission l c (0e) by the to- 
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FIG. 4: (color online) Angular distribution of the R6G emis- 
sion from a polystyrene opal with a lattice parameter a — 178 
nm. The black pentagons indicate the intensity measured at 
the spectral maximum (ui = 17860 cm -1 or A = 560 nm), 
the blue curve is the calculated escape function with Fresnel 
internal-reflection coefficient. Measured intensities corrected 
for the detection efficiency of the set-up are displayed as red 
stars. All data are normalized at 8 e = 0°. 



tal emission spectrum I to t and correcting for the angular 
aperture of the collecting lens Ql- For the frequency u 
= 15000 cm -1 , below the stop band, the distribution fol- 
lows the Lambertian distribution and is similar to the 
exit distribution from the nonphotonic sample (Fig. 0}. 
For the frequencies above the red edge of the stop band 
we observe strongly non-Lambertian behaviour. For the 
frequency u> = 16500 cm -1 emission is suppressed rela- 
tive to the Lambertian distribution in the range of the 
exit angles from 8 e = 0° to 20°. This range moves to 
larger exit angles for the frequency ui — 17400 cm -1 in 
qualitative agreement with Bragg's law. For u> = 18350 
cm -1 the suppression observed around 8 e = 40° is pre- 
ceded by a considerable increase of emission in the angle 
range 0° to 20°. A qualitative explanation of this effect is 
as follows. Some escape directions are blocked by inter- 
nal Bragg diffraction, and diffusion eventually distributes 
this back-reflected light along the remaining directions. 
Thus, light is more likely to exit the crystal along these 
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FIG. 5: (color online) R6G emission from polystyrene opal 
with lattice parameter a = 365 nm. The scatter plots repre- 
sent the measured spectra corrected for the angular aperture 
of the collecting lens Ql and divided by the total emission 
spectrum Itot(u). The calculated escape functions are plot- 
ted with solid curves. The black solid curve and squares are 
for 9 e = 0°, red dashed curve and triangles are for 8 e = 30°, 
green inverted triangles and solid curve are for 6 e = 45° and 
blue circles with dash-dotted curve are for 8 e = 60°. 
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FIG. 6: (color online) Intensity ratios I c (6 e )/ hot corrected 
for the lens aperture Ql as a function of the exit angle t9 e for 
a polystyrene opal with a lattice parameter a = 365 nm for 
frequencies lo = 15000, 16500, 17400 and 18350 cm -1 (blue 
inverted triangles, black squares, red circles and green tri- 
angles). The corresponding curves represent the calculated 
escape distributions (for reduced frequencies uia/2nc = 0.55, 
0.6, 0.64 and 0.67, respectively). 

allowed directions. From Figures |3] and we conclude 
that the escape function is in excellent agreement with 
the measured angle-dependent spectra. To the best of 
our knowledge, the current work provides the first quan- 
titative modelling of spontaneous emission spectra in 3D 
photonic crystals. 



B. Spontaneous emission from quantum dots in 
titania inverse opals 

Titania inverse opals possess a larger relative width 
of the first-order stop band (Aw/u « 0.16) than the 
polystyrene opals owing to their inverse structure and 
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high refractive index contrast (m = 2.7). The concomi- 
tant large modifications of the LDOS makes the inverse 
opals very attractive for control of propagation and spon- 
taneous emission of light 0,^3- Figure [7| shows emission 
spectra of CdSe quantum dots in a titania inverse opal 
with lattice parameter a = 500 nm for selected detec- 
tion angles 9 e . No significant changes in the spectral 
shapes due to internal Bragg diffraction are observed, 
because the relative spectral width of the light sources 
(Aw/w < 0.06) is considerably smaller than the width 
of the stop band of the photonic crystal. This shows 
that the escape distribution P(u,9 e ) does not vary sig- 
nificantly within the frequency range of the quantum dot 
spectrum. In contrast, a strong angular dependence of 
the emission intensity is apparent in Fig. [7] As a conse- 
quence, effects of Bragg diffraction are most convincingly 
observed by recording the angular dependencies at the 
spectral maxima of the emission spectra. 

In Figure Ufa) we present escape distributions from 
titania inverse opals with lattice parameters a = 370 
and 420 nm at frequency to = 16390 cm -1 , and with 
a lattice parameter a — 500 nm at u = 15870 cm -1 . 
Both measured (symbols) and calculated (curves) values 
are shown. For the crystal with the lattice parameter 
a = 370 nm, for which the center frequency of the quan- 
tum dot spectrum lies below the stop band, the escape 
function follows the Lambertian distribution. A large 
deviation from the Lambertian distribution is observed 
for the quantum dot emission from the crystals with the 
other two lattice parameters. In the crystals with a = 420 
nm, the emission is strongly reduced in the range of the 
angles from 8 e — 0° to 35°, and it is enhanced at higher 
exit angles. For the crystals with a = 500 nm, the sup- 
pression is shifted to the range of 9 e — 20° to 45°, as 
expected for photonic gaps at higher reduced frequency 
(a/A = Lua/2nc), and in excellent agreement with our 
theoretical predictions. The stop-band ranges are no- 
ticeably wider than that in the case of the polystyrene 
opals (Fig. EI), which is due to a wider frequency range of 
the stop band in the titania inverse opals. To the best of 
our knowledge, this is the first demonstration of photonic 
crystal bands in the emission spectra of confined excitons 
in quantum dots. 

Figure IBfb) shows the photonic band structure for a 
titania inverse opal. The hatched regions indicate the 
stop band caused by Bragg diffraction by (111) lattice 
planes. In the angular range from 9 e — 35° to 55°, multi- 
ple Bragg wave co uplin g from (111) and (200) diffracted 
waves takes place [la . l36j . The horizontal bars repre- 
sent the reduced center frequencies of the quantum dot 
emission from the inverse opals with the lattice param- 
eters a = 370, 420, 500, 580, and 650 nm. The colours 
of the bars indicate the measured values of the escape 
function P{u>,9 e ). For reduced frequencies around the 
stop bands, it is seen that inhibited escape probability 
appears in the angular ranges of the stop bands, whereas 
enhanced escape is found outside the stop bands. Hence, 
the photonic crystals are seen to 'funnel' light along cer- 
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FIG. 7: (color online) Emission spectra of CdSe quantum dots 
in a titania inverse opal with a lattice parameter a = 500 nm. 
The black curve is obtained at 6 e = 0°, the red dashed curve 
at e — 15°, the green dash-dotted curve at 8 e = 25°, and the 
blue short-dashed curve at 9 e = 50°. 



tain allowed directions. 

For experiments on quantum dot emission in photonic 
crystals, the good agreement between the experimen- 
tally obtained escape distributions and the calculated 
ones confirms that the light emanating from inside the 
crystals is diffuse. It also confirms that the observed 
emission is dominated by sources inside the bulk of the 
crystal. We can exclude that light sources on the sam- 
ple surface contribute significantly: their emission would 
give rise to an angle-independent component of the in- 
tensity that is not observed. Furthermore, observation of 
stop bands in emission spectra is important for successful 
lifetime experiments or other quantum-optical studies of 
light sources in photonic crystals. The stop bands are ev- 
idence that the emission from the light sources is strongly 
coupled to the photonic crystals, and are a prerequisite 
for time-resolved experiments of changes in the emission 
decay rate caused by a modified LDOS |Toj . 

Based on the close accordance of the experiments with 
the model, we can extract the angle-averaged internal- 
reflection coefficient R(u>). Figure [5] shows that R(u>) is 
as large as 50% for the titania inverse opals and up to 
20% for the polystyrene opals. The internal-reflection 
coefficient varies strongly with frequency in contrast to 
the frequency-independent I? in random media. The co- 
efficient increases with solid angle for Bragg reflection, 
starting from the low-frequency edge of the L-gap. For 
the opals, the maximum R occurs at the high-frequency 
edge of the L-gap, where the reflecting stop bands ex- 
tend over the largest solid angle j34[. For the inverse 
opals, the maximum R occurs at higher reduced fre- 
quencies in the rang e of multiple Bragg wave coupling 
(LoajlTtc ~ 0.85) [15|. The shoulder near 1.0 is attributed 
to the inclusion of (200) reflection condition in our model. 
In a more elaborate escape-model with additional diffrac- 
tion conditions, we may expect additional peaks in the 
angle-averaged reflection coefficient at even higher fre- 



quencies. Since the inverse opals interact stronger with 
light than the opals, their stop bands are wider and hence 
the angle averaged reflectivity is larger, in agreement 
with our observations. 

A closer consideration of the reflectivity coefficients 
can serve to optimize the spontaneous emission yield of 
light sources (atoms, dyes or quantum dots) embedded 
inside thick photonic crystals (L > I). Such an opti- 
mization can be achieved either via the excitation of the 
sources, via their emission, or both. First, the excitation 
efficiency can be increased by realizing that increased es- 
cape probability also implies an increased probability for 
excitation light to enter a photonic crystal. Thus, by 
tuning an excitation beam to frequencies and angles of 
high escape, the combined action of diffusion and Bragg 
diffraction retain relatively more excitation light inside 
the sample, thus increasing the probability for sponta- 
neous emission of the embedded light sources. Secondly, 
spontaneously emitted radiation is efficiently channelled 
out of the sample along particular directions. This oc- 
curs when the lattice parameter of the photonic crystal 
is chosen such that the emission frequencies are in the 
range of enhanced R(ui). A clear example of enhanced 
escape is apparent in Fig. at 9 e = 0° near 18000 cm -1 . 
In the ultimate case of a photonic band gap, it has even 
been predicted that the diffuse emission is extremely di- 
rectional, see Ref. [2l|. Thirdly, one can envision situa- 
tions where both excitation and emission are enhanced: 
In Fig. Efb) enhanced escape probability occurs both at 
Loa/2nc ~ 0.8 and 9 e — 0° and at toajl-Kc ~ 0.7 and 
6 e = 50°. Thus, by tuning the excitation to the former 
condition and the emission to the latter, the spontaneous 
emission yield is expected to be enhanced by at least a 
factor of two. Further improvements should be feasible 
in photonic crystals with even larger R. 



V. CONCLUSIONS 

We have presented experimental data on angular re- 
solved emission from light sources embedded in efficient 
3D photonic crystals. The experiments were compared 
in detail to a recently developed model of light transport 
in real photonic crystals that are influenced by disor- 
der. Our model is based on diffusion of light due to scat- 
tering (disorder) combined with angle- and frequency- 
dependent internal reflections (order). Good quantita- 
tive agreement between experiment and theory confirms 
that the details of the emission spectra are determined 
by the intricate interplay of order and disorder. Prop- 
erties of the stop bands, such as their frequency range, 
magnitude, and angular dependence, are extracted from 
the experiment by analyzing the emission escape func- 
tion. The enhanced escape probability for emission along 
directions outside the stop bands is explained by the 
angle-averaged internal-reflection coefficient R(uS). The 
diffuse and angular-dependent nature of light escaping 
from the photonic crystals proves that the light comes 
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FIG. 8: (color) (a) Escape distribution as a function of the exit angle 8 e for titania inverse opals with lattice parameters a — 
370 and 420 nm at u> = 16390 cm -1 (A = 610 nm, black squares and red circles, respectively), and with a lattice parameter a 
= 500 nm at u — 15870 cm -1 (A = 630 nm, green triangles). The corresponding curves represent the calculated distributions, 
(b) Photonic band structure for the inverse opals (black dotted curves). The hatched regions indicate the stop band caused 
by Bragg diffraction by (111) lattice planes. The horizontal bars represent the reduced center frequencies of the quantum-dot 
emission from the crystals with lattice parameters (bottom to top) a = 370, 420, 500, 580, and 650 nm. The colours of the bars 
indicate the values of the escape function P(oj,8 e ) obtained from the measurements at exit angles shown by the black stars. 



from emitters inside the crystals. By measuring the es- 
cape functions of the quantum dot emission from the ti- 
tania inverse opals, we have for the first time revealed 
clear stop bands in the quantum dot emission spectra, 
confirming that the confined excitons experience optical 
confinement. The quantitative agreement between ex- 
periment and theory demonstrates that light propagation 
and spontaneous emission in real 3D photonic crystals is 
well understood. 
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VII. APPENDIX: MODEL OF DETECTION 
EFFICIENCY OF EMISSION SET-UP 

The aim of this appendix is to explain the difference 
between the calculated escape distribution and the mea- 
sured intensity illustrated in Fig. This difference ap- 
pears because the width of the projection of the collec- 
tion optics (the spectrometer slit) on the sample increases 
with the angle 9 e . The sample surface is placed in the 
focus of the collecting lens L\, the spectrometer slit is 
in the focus of the imaging lens L 2 (cf. Fig. P). The 
only emission collected emanates from the region con- 
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FIG. 9: (color online) The angle-averaged internal reflectiv- 
ity R(u>) for the polystyrene opals (black curve) and for the 
titania inverse opals (red dashed curve) according to the dif- 
fusion model, in which only the first-order Bragg diffraction 
is taken into account. R(io) determines the enhancement of 
the escape probability outside a stop-band direction. 



Emission spot 
Slit projection 




FIG. 10: (color online) Real-space cartoon of the projection of 
the spectrometer slit (striped rectangle) on the sample surface 
overlapped with the emission spot (in blue). The quarter of 
the slit projection (red dashed rectangle) has one of the cor- 
ners in the emission-spot center. The width of the projection 
(along x) increases at larger detection angles 8 e . 
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fined by the slit projection on the sample surface, see 
Fig. This means that the spectrometer collects light 
from a larger region on the surface at larger detection 
angles, and that the measured angle-dependent intensity 
should be corrected for the detection efficiency of the set- 
up. The detection efficiency is modelled as a ratio D(fi e ) 
of the intensity B(fi e ) collected from the surface region 
within the slit projection (Fig. I1U[1 at detection angle 
9 e = cos _1 (^ e ) to the intensity B(/i e = 1) collected at 
normal angle: 

D ^ = Wfj, = i) ' B ^ = J J I ( x >y) dxd y- 

(7) 



We take into account that the integration runs over a 
quarter of the slit, as x(fi e ) = Xo/(A e is the half-width of 
the slit projection, :ro and yo arc the half-width and the 
half- height of the slit projection at fi e — cos(9 e ) — 1, re- 
spectively. Typical values of xq and yo in the experiments 
are 50 fim and 1 mm. It is assumed that the intensity 
of diffuse light on the sample surface around the pump 
beam varies as J(r) oc l 2 /(l 2 + r 2 ), where r 2 = x 2 + y 2 is 
the distance from the axis of the pump beam along the 
sample surface, and I is the mean free path of light in the 
sample. 
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